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Abstract Magnetic field diagnostics of the transition region from the chromo- 
sphere to the corona faces us with the problem that one has to apply extreme 
UV spectro-polarimetry. While for coronal diagnostic techniques already exist 
through infrared coronagraphy above the limb and radio observations on the 
disk, for the transition region one has to investigate extreme UV observations. 
However, so far the success of such observations has been limited, but there are 
various projects to get spectro-polarimetric data in the extreme UV in the near 
future. Therefore it is timely to study the polarimctric signals wc can expect for 
such observations through realistic forward modeling. 

We employ a 3D MHD forward model of the solar corona and synthesize the 
Stokes I and Stokes V profiles of Civ (1548 A). A signal well above 0.001 in 
Stokes V can be expected, even when integrating for several minutes in order 
to reach the required signal-to-noise ratio, despite the fact that the intensity in 
the model is rapidly changing (just as in observations). Often this variability of 
the intensity is used as an argument against transition region magnetic diag- 
nostics which requires exposure times of minutes. However, the magnetic field 
is evolving much slower than the intensity, and thus when integrating in time 
the degree of (circular) polarization remains rather constant. Our study shows 
the feasibility to measure the transition region magnetic field, if a polarimctric 
accuracy on the order of 0.001 can be reached, which we can expect from planned 
instrumentation. 

Keywords: Sun: transition region — Sun: corona — Sun: UV radiation - 
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1. Introduction 

There is a general consensus that the outer atmosphere of the Sun (and othcr 
cool stars) is heated by one or s everal mechanisms related to the magnetic field 
(e.g. Schriiver and Zwaan . 2000l) . Despite the pivotal importance of the magnetic 
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field for our understanding of the corona, actual measurements of the magnetic 
field are scarce — mostly we have to rely on ext rapolations of the magnetic field 
from the photosphere fe.o. jDe Rosa etcd. . 20091 ). Because the extrapolations are 
based on assumptions that might or might not be fulfilled in regions of interest 
in the corona, there is the need to actually measure the coronal magnetic field. 

Some measurements in the corona have been performed in active region loops 
above the limb using the Z eema n effect for infrared coronagraphic observations 



( Lin. Penn. and Tomczvk , 2000() . These, as well as radio measurements (e.g. 



White! . 2005 ). suffer from low spatial and temporal resolution. A very promising 



project for diagnostics of the magnetic field in the upper chromosphere is the 
rocket experiment chr omospheric Lyman-alpha spectro-polarimeter (CLASP; 



Kobavashi et al. . 2012[) . planned to be flo wn in Dec 2014. This is based on di 



agno stics using the Hanle effect wiin Ly-a ( Truiillo Bueno. Steoan. and Casinil . 



20111) . 



The first attempt to measure the magnetic field on the solar disk at high 
resolution in the transition region from the chromosphere to the corona was 
done using the Sola r Maximum Mis s ion ( SMM) ultraviolet spectrometer and 
polarimeter (UVSP; fWoodgate et al. . 1980l ). Reaching a polarimetric accuracy 



of jus t below 1% in the ClV line at 1548 A ( Henze et all Il982t iHagvard et al 



1983 ) gave no conclusive results, except maybe in sunspots ( LitesL 2001 ). As it 



will become clear in this manuscript, an accuracy of 0.1% would be needed to 
get useful information on the magnetic field with the C IV line. This accuracy is 



provi ded by the solar ultraviolet magnetograph investigation (SUMI: fWest et al 



20041 ) that has been flown on a rocket twice. The C iv data from the latest flight 



in summer 2012 await calibration and analysis. 

The present study will present a forward model which provides synthesized 
polarimetric data as they are recorded, e.g., by SUMI. This will allow a di- 
rect comparison between model and observations and will (hopefully) provide 
some guidance for the interpretation of the acquired polarimetric observations. 
A 3D model is prerequisite to investiga te the transition r egion, because of its 
highly complex spatial structuring (e.g. IPeter . |2000L l200ll) . Consequently, this 
study is based on a 3D magnctohydrodynamic (MHD) numerical experiment 
that provides temperature, density, velocity and of course the magnetic field in 
the cor ona above a small active r egion . This m odel produces a loop-dominated 
corona (|Gudiksen and Nordlundl . l2002i |2005allbh. In a statistical sense it repro- 
duces various observational properties ( Peter. Gudiksen. and Nordlundl. 2004 . 



2006 ) , in particul ar the persistent transition region redshifts (jPeter and Judge 



1999t ; |Peterlll999h . Based on the success of these models the present study goes 



one step further to investigate not only the profiles of the emission lines, but also 
the circular polarization due to the Zeeman effect. Another forward modeling 
investigation of the polarization of the emergent radiation in an MHD model of 
the extended solar atmo sphere, but fo r the ( optically thick) hydrogen Ly-a line, 



has been carried out bv lStepan et al\ (|2012f ) paying particular attention to the 



linear polarization produced by scattering processes and the Hanle effect. 

Once the instruments provide the required polarimetric sensitivity, we will 
have to perform a proper interpretation of the data. The extreme UV lines 
formed in the transition region (or in hotter parts of the atmosphere) originate 
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from a spatially complex volume, comparable to highly corrugated surfaces. This 
is a fundamental difference to photospheric magnetic field observations, where 
the height of the source surface remains at a roughly constant altitude, within 
one barometric scale height. This more complex source region of transition region 
lines complicates possible inversions enormously. Furthermore, the high temporal 
variability of the intensity will be a major problem for the interpretation of the 
polarimctric data. The forward-model approach presented in this study pro- 
vides insight into these problems by accounting for all the spatial and temporal 
complexity. 

In Section^] the basic concept of the 3D MHD model and the spectral syn- 
thesis including the calculation of the Stokes V profile will be introduced. Based 
on this in Section[3] we will present how to synthesize the observable quanti- 
ties and present some sample Stokes V profiles in Scction|4l The observational 
requirements for the exposure times will be discussed in Section[5] along with 
the spatio-temporal variability found in observations and in synthesized model 
data. In Section^ we will construct a realistic Stokes V observation of C IV and 
show that we can afford comparably long exposure times for the observations — 
and why. Finally, in Section[7]we will discuss some simple inversions and their 
reliability, before we conclude the paper with Section[Sl 



2. Synthetic spectra from 3D MHD model 



In general, the state of polarization of the light can be described by the Stokes 
vector (/, Q, U, V). Stokes I represents the integral over all polarization states, 
while Stokes V is the difference of right- and lcft-circularly polarized light, and 
hence carries information on the longitudinal component of the magnetic field 
through the Zccman effect. For an observation near disk center and assuming 
that the magnetic field will be predominantly vertical, we can expect Stokes Q 
and U, which characterize the linear polarization, to be much weaker than Stokes 
V, just as found in photospheric observations in the visible. Because it will turn 
out that already the Stokes V signal will be at the edge of observability, this 
study will concentrate on Stokes / and V only. 



2.1. 3D MHD model 



In order to calculate the Stokes profiles, one needs the temperature T, density n, 
velocity vector v, and magnetic field vector B. These are provided by a 3D MHD 
model that solves for the induction equation, the conservation of mass, and the 
momentum and energy balance. Most importantly, the energy balance has to 
include heat conduction and radiative losses. This is pivotal to set the proper 
coronal pressure and therefore a prerequisite to synthesize coronal emission lines 
which are very sensitive to the temperature and density. 

The MHD model u s ed for this study has been published already by 
Gudiksen and Nordlundl (l2002Ll2005allbh . and an analysis of this model m terms 



of Doppler sh ifts and emission measure and a comparison to observations was 
presented by IPeter. Gudiksen. and Nordlundl ( 20041 l2006h . In the model the 
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plasma is heated through Ohmic dissipation of currents that are induced by the 
braiding of magnetic field lines through the horizontal photospheric granular 
motions. The good match to the observations showed that this model provides a 
realistic way to describe the corona in an active region, accounting for the spatial 
and temporal variability. Since then further models of similar type solidifie d 
these results discussing details of the heat input ([Bingert and Peterl . l201lll2013l) . 
providing further insight into the persistent transition re gion Doppler shifts 
( Hansteen et al. 1 2010t Zacharias. Bingert. and Peter , 2011b ) , transient events in 
the corona ( Zacharias. Bingert. and Peteru2011a ). or the constant cross section 



of loops ( Peter and Bingert! 2012t) All these results give good confidence that 



this type of model will also allow for a reliable and realistic determination of the 
Stokes V profiles in a transition region extreme UV line. 

2.2. Intensity spectra: Stokes / profiles 



To calculate the Stokes profi l e 7(A) we follow exactly the procedure of 
Peter. Gudiksen. and Nordlundl ( 2004L [2006). Under the assumption of ioniza- 
tion equilibrium we calculate the emissivity (energy loss per time and volume) 
at each grid point in the compu tational domain in the Civ (1548 A) line using 
the Chianti atomic data base ( Dere et all . Il997t lYoung et all . l2003h . In order 
to avoid aliasing effects we interpolate the mesh in the vertical direction. The 
line width is assumed to be the thermal width, A Ad, and the Doppler shift is 
given by the line-of-sight component of the velocity vector (here the vertical 
component). This provides a line profile I(X) at each grid point. 



2.3. Stokes V profiles: weak-field approximation 

In a magnetic field a spectral line will be affected by the Zeeman effect. The 
Zeeman splitting is given through 



AAr = 



Air cm. 



9$B 



los ; 



(1) 



with the elementary charge e, the speed of light c, the electron mass m e , the rest 
wavelength Ao, the effective Lande factor g and the component of the magnetic 
field along the lone of sight B\ os . 

In the weak-field limit the Stokes V signal is given through thr ough a Fourier 
expansion up to first order of the Stokes I line profile (e.g., Stenflol . 1994 . 
Section 11.9), 



V(X) 



-AA, 



dl_ 



4.67-10" 13 5(A [A]) 2 B los [G] 



dl 

a a [A] ' 



(2) 



In this study we will concentrate on the C iv line at 1548 A when investigating 
the transition region magnetic field, as it has a high diagnostic potential: it is a 
strong line at a long wavelength Ao (for transition region extreme UV lines) with 
a decent effective Lande factor g. This combination provides the best potential 
among the transition region lines. 
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The effective Lande factor for the 1548 A lin e ( 2 P 3/ o -±fSi/ 2 ) is 5=7/6^1.167 
(for details on the calculation of g see e.g. IStenflol . 11994 Sects. 6.4 & 6.5). 



While the 1550 A line of the Civ doublet has a slightly higher effective Lande 
factor (5=4/3), its radiance is lower by a factor of about two, which is why we 
concentrate on the 1548 A line. 

Even for a very strong magnetic field in the transition region of 1000 G the 
splitting for the Civ lines would be below AAs<0.15pm corresponding to less 
than 0.3 kms -1 in Dopplcr shift units. This can be considered as an upper limit. 
Consequently, in the case of C iv the Zeeman splitting is much smaller than 
the Doppler broadening of the line, AAb^cAAd, the latter being about 6 pm 
corresponding to 12kms _1 at the line formation temperature of Civ of about 
10 5 K. Thus the application of the weak-field limit for Equation ([2]) is justified. 

From Equation we can roughly estimate the Stokes V/I signal to be ex- 
pected in C iv. For this we assume a magnetic field of 100 G (above a pore or a 
strong network patch) in the transition region some 3 Mm above the photosphere. 

Then, for a line width of the order of the thermal width the expected signal is 
about V/I f» 0.002, which would be measurable with planned instrumentation 
(cf. SectioniJ). 

With the line-of-sight magnetic field from the 3D model and the Stokes 
I profile from Section l2~2l we can compute the Stokes V profile according to 
Equation ([2]) at each grid point on the interpolated mesh of the computational 
domain. 

In this study we will concentrate on the Civ line formed at about 10 5 K. It 
could be easily also done for any other extreme UV emission line, e.g., the Ovi 
line at 1032 A formed at some 300 000 K, or the Mgx line at 625 A formed at 
10 6 K to study the coronal magnetic field. However, at the shorter wavelengths 
the Stokes V signal will be weaker. Furthermore, at higher temperature, on 
average found at higher heights, the magnetic field will be weaker. Thus to 
detect a Stokes V signal from these lines the sensitivity of future instruments 
would have to be well below 0.1%. 



3. Synthetic observation and diagnostics 

For this study we will restrict the discussion to a part of the computational 
domain covering about 27 x 27 Mm 2 in the horizontal directions. This region of 
interest contains a magnetic structure at the bottom boundary that resembles 
a pore on the real Sun. The (vertical component of) the magnetic field at the 
bottom boundary, i.e. the photosphere, in this area is shown in Figurc[T](a) . In 
this region the photospheric magnetic field is predominantly of one sign. 

We will investigate synthetic extreme UV observations for which the line of 
sight is aligned with the vertical. For this we integrate the Stokes /(A) and V(A) 
profiles along the vertical direction, 

7 los (A) = / /(A) dz ; Mo S (A) = / V(X) dz . (3) 
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Figure 1. Snapshot of a view when looking from straight above at the computational box 
(along the vertical direction). Panel (a) shows the line-of-sight magnetic field at the bottom of 
the box defined through the boundary condition. Panel (b) displays the synthesized emission 
of the C IV line, integrated along the line of sight. The line-of-sight component of the transition 
region magnetic field at the height of the formation of C IV is seen in panel (c). Panel (d) shows 
the resulting integrated Stokes V signal normalized by the line intensity (see Eq.|4]l. The field 
of view shown here is about 27 X 27 Mm 2 and represents part of the whole computational 
domain. The two yellow crosses denote the locations for which the sample Stokes V spectra 
arc plotted in FigurcT3h. The red and green boxes show the regions over which realistic Stokes 
profiles are synthesized (Figure|8jl and in which the temporal variation is analyzed (Figure|6j|. 
The long blue bar indicates the line along which the spatial variation is analyzed (Figurc[7]). 



This will produce observables as they would be found in actual observations of 
the Sun close to disk center when employing an extreme UV spectro-polarimeter. 
These data can then be analyzed in the same way as actual observations of the 
Sun would be handled. For instance, one can derive maps of the total line inten- 
sity, Doppler shifts, or the wavelength-integrated Stokes V signal. In particular, 
we define the total intensity and the integrated Stokes V signal as 

T = jl los (\)d\ ; overlineV = ^ J \V los \(X) dX . (4) 

The latter gives the fraction of the (unsigned) Stokes signal along the line of 
sight when compared to the total line emission. In Figures[T|D andQJl we show 
the resulting maps of / and V for the part of the computational domain under 
investigation here. 

Having access to the three-dimensional distribution of the C iv emission in the 
computational domain, at each horizontal location (x, y) we can derive at which 
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(a) C IV (1548A) line formation height: 3... 11 Mm 




Figure 2. Line-formation height of Civ. The surface in panel (a) shows the height where 
C IV is originating. Displayed on that surface is the emission from C IV to show the (inverse) 
relation of emission and formation height. The horizontal extent of the field of view is the same 
as in Figure^ Panel (b) presents the distribution of formation heights in those areas where 
the Stokes signal V as defined in Equation Q exceeds 10 — ' ? . 



height z the main contribution of C iv is originating from. This defines the line- 
formation height, which is visualized in Figure^. As noted before, this is a highly 
corrugated surface. Typically, the line-formation height is the lowest where the 
intensity is the highest - where the heating is particularly high the transition will 
move to lower heights and produce more emission. Even when considering only 
regions with considerable field strength, where the Stokes signal V as defined in 
Equation (QJ exceeds 10 -3 , there is a wide range of heights of the source region 
of C IV. The histogram in Figure^}) shows that there the formation of C IV is 
mostly found between 3 Mm and 4 Mm, sometimes reaching up to 6 Mm. This 
corresponds to several chromospheric pressure scales height («300km). 

Having the line- formation- height, we can now extract the magnetic field at 
the source region of C iv. We plot the vertical component of this, i.e., the line-of- 
sight component, in FigurejlJ;. This map of the transition region magnetic field 
is not just a horizontal cut, but shows the magnetic field at the actual height 
of the main transition region emission for each horizontal location, i.e., on the 
corrugated surface shown in Figure^- 

This surface of the maximum contribution of C iv coincides with the lo- 
cation where the temperature jump into the corona is found. However, the 
spatial structure of the transition region is even more complicated than this. 
The emission of the optically thin C iv line is not restricted to this surface, 
but pockets of 10 5 K cool plasma are also found higher up in the (generally) 
hotter volume of the corona, as can be seen from the images and movies of 
Peter. Gudiksen. and Nordlund (l2006h . One should keep this in mind when in- 



terpreting the data - along the line-of-sight various structure can contribute to 
the emission seen in transition region lines. 

Comparing the photosphcric and transition region magnetic field in Figures^ 
and[Tl:, it is clear that the transition region field is much smoother and no longer 
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wavelength shift [km/s] Stokes V signal: JlV* los ldX / jl lot dX 

Figure 3. Synthesized Stokes V profiles and histogram of the distribution of Stokes V signals. 
Panel (a) shows two sample Stokes V profiles normalized to the peak line intensity for the loca- 
tions marked by crosses in Figure^ (the stronger signal shown as a solid line corresponds to the 
right cross). Panel (b) displays the distribution of Stokes V signals as defined in Equation Q 
for the field of view shown in Figure^ About one third of the area shows a Stokes V signal 
above 0.1 %. Alternatively one could also plot the histogram of the peak values of the Stokes 
V profiles, which would show a very similar result. Sec SectionfJ] 



shows any signs of (small-scale) mixed polarities. The opposite polarities in the 
photosphere are a couple of Mm apart, so at a height of 3 Mm and higher, 
where C iv is formed, all these mixed polarities are closed already. Similar to a 
simple potential field expansion, the 3D MHD model shows the expansion and 
smoothing of the magnetic field with height. 



4. Stokes V line profiles 



Two examples of Stokes V profiles are shown in Figure^. The sample from the 
location close to the center of the magnetic field concentration in the photosphere 
shows a peak-to-peak signal in Stokes V of C iv normalized to the peak intensity 
of almost 1%. This profile is among the strongest Stokes V signals in the ficld- 
of-vicw (right cross in Figure^. The signal is symmetric and close to a Stokes 
V signal in an idealized situation. In regions where the magnetic field is weak, 
but where still significant emission in ClV is present (left cross in Figured]), the 
Stokes V signal typically is weaker than 0.1%. 

The distribution of the C iv Stokes V signals in the field-of-view displayed in 
Figure[3p reveals that about one third of the are a shows a signal ab ove 0.1%. This 



is within the detection limit of the SUMI rocket (TWest et all 120041) . Likewise, the 



future rocket mission CLASP (jKobavashi et all l2012|) will achieve this accuracy 



in the extreme UV (albeit at Ly-a) . Furthermore the space-based observatory 



SolmeX that was proposed to ESA (IPeter et all 120121 ) would have included a 



spectro-polarimeter with the required capabilities for C iv observations (see also 
Section©. 

The majority of the area in the field-of-view visible in FigureQ] shows signals 
much smaller than 0.1% (cf. Figure^) that will not be detectable. Mostly these 
weak signals originate from regions with very low emission in C iv. So besides 
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5 10 

observation time [min] 



Figure 4. Noise-to-signal ratio for a photon-counting detector (SUMER/SOHO) from an 
actual observation. The data show Si IV at 140.3 nm observed in the quiet Sun network. After 
12 minutes 10 6 counts were reached and the noise level dropped to 10 -3 . See Section fe-ll 



the weak Stokes V signal also the low intensity would prevent the detection of 
a signal here because of the limitation in signal-to-noise. 

The distribution in Figured shows the Stokes signal as defined in Equa- 
tion (QJ, i-e., the total unsigned Stokes V normalized by the total intensity. The 
distribution for the peak values of Stokes V normalized by the peak intensity 
would show a very similar distribution. 



5. Variability of the transition region intensity and magnetic field 
5.1. Observational requirements for exposure times 

Based on the discussion of Figures[T}l and [3]d in Section|4l and considering the 
order-of-magnitude arguments in Section [231 it is clear that instruments would 
have to detect a Stokes V signal (normalized to I) of 0.001 or better. This 
implies that the signal-to-noise ratio of the recorded data has to better than 
1000. Assuming Poisson statistics this leads to the conclusion that the detector 
has to acquire at least 10 6 counts (or 10 7 counts for a 3cr detection). This will 
set a limit of the required exposure time. 

As a simple experiment one can investigate a time series of actual obser- 
vations to get an estimate of the exposure time. The most recent instrumen t 
that recorded high-quality solar data in C iv is SUMER dWilhelm et al.l . ll995h . 



Because of the lack of a well suited time series in C IV here data from Si IV are 
used, a line that forms at similar temperatures as C IV. The two lines share the 
major observational properties, e.g., variability, line shifts, contrast, etc. 

To investigate the noise level, a time series with high temporal cadence is 
analyzed and the acquired counts are accumulated. This is possible with a 
photon counting device as SUMER. Assuming Poisson statistics, which is a good 
approximation for the SUMER detector, the noise-to-signal ratio is derived. This 
is simply the inverse of the square root of the accumulated counts. In FigureH] 
this noise-to-signal ratio based on the accumulated counts is shown as a function 
of the observational time. As expected the noise level generally drops in time. 
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distance along the slit [Mm] time [min] 



Figure 5. Observed actual spatial and temporal variability of the transition region. Example 
of SUMER data showing the intensity of the Ovi line at 1038 A observed on 16 May 1998. 
They represent a region with strong network activity. The left panel (a) displays the spatial 
variation (along the slit), with a pixel size of 1" (crosses), i.e. at the highest possible spatial 
sampling, at an exposure time of 10 s. The right panel (b) shows the temporal variation in a 
single l"xl" spatial pixel (crosses). The bars in the panels indicate the average uncertainty 
for the line intensity (Poisson noise). The crosses indicate the spatial pixels (left) and times of 
exposures (right). The temporal and the spatial variability is down to the resolution limit of 
SUMER. See SectionlOl 

The details depend on the time-variable emission, of course. For the particular 
example shown in Figure^] after some 12 minutes a noise level of 0.001 is reached. 

Therefore, to reach a signal-to-noise level to detect a Stokes V signal of 0.001 
one would have to integrate for some 12 minutes, too. The data shown here are 
for a bright patch in the network, so for regions with higher magnetic field, e.g., 
in the vicinity of a pore, we can expect higher emission in the transition region 
lines, and thus a shorter exposure time would be sufficient. Also a more modern 
instrument could have a higher throughput and larger aperture, which would in 
part be counterweight by the additional optics needed to analyze the state of 
polarization of the incoming light. So as a very rough estimate one might need 
exposure times of the order of minutes. 

In the light of the high variability of the transition region emission this 
poses the question if Stokes V observations with such comparably long exposure 
times are meaningful. After presenting the transition region variability in actual 
observations (Section IO)) and in the synthesized model data f Section lOl) the 
discussion will turn to this question in Scction[B] 

5.2. Actually observed transition region variability 

Just as an example, Figure[5] shows the temporal and spatial variability of the 
transition region emission. The spectra were acquired in the quiet Sun in a 
region with strong network activity. The spatial variation (panel a) shows a 
35 Mm («48") long cut through two network patches with enhanced emission. 
Besides this spatial variation on a scale of some 10 Mm, also variations down to 
the resolution limit (1" spatial sampling) can be seen. It has to be emphasized 
that this pixel-to-pixel variation is not noise, but represents real variability not 
fully resolved by the instrument. 
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10 12 14 
time [min 

Figure 6. Temporal variability of synthesized transition region emission and magnetic field. 
The top panels (al) and (bl) show the temporal variation in the two 1.6 X 1.6 Mm 2 squares 
labeled in FigurcTlk as (a) and (b). The thick solid lines show the variation of the average 
intensity in the respective square, while the thin lines with the crosses show the variation at 
the central grid point in that square. The bottom panels (a2) and (b2) show the variation of 
the magnetic field at the line formation height of C IV, i.e., the source region of the transition 
region emission. Again, this is done for the average within the respective square (thick) and 
the central grid point (thin with crosses). The dashed lines indicate the magnetic field inverted 
from the average Stokes V signal as discussed in Section]?] See Scction[fO] 



For the temporal variation (Figure^) one can identify fluctuations also down 
to the sampling of the instrument (here 10 s). Again, these fast fluctuations are 
real, and consistent with the very short cooling times in the transition region. 
Stronger fluctuations well above a factor of two occur in time scales of minutes 
that are omnipresent on the Sun. In this example two cases, one shorter and 
one longer can be identified. Thes e brightenings can be classified as b l inkers 
that have been studied abundantly (Harrison et al. , 19991 Harrison et al , 2003; 



Peter and Brkovid 120031: iBrkovic and Peterl . 1200 J ) 



This strong variability in time (and space) raises the question to what ex- 
tent spectro-polarimetric data will be useful to investigate the transition region 
magnetic field with instruments that will have limited spatial and temporal 
resolution. This can be explored by observations synthesized from a model, where 
the magnetic field is known. 

5.3. Spatio-temporal variability in synthesized observations 

Examples for the temporal variability of the transition region emission syn- 
thesized from the 3D MHD model are shown in Figure^] (top panels) for the 
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squares labeled (a) and (b) in Figure[T] These regions are selected to be in regions 
near the main polarity with medium (a) and high (b) intensity in C iv. In the 
figure the variation integrated over an area of 1.6x1.6 Mm 2 (roughly 2"x2") is 
plotted along with the intensity of the simulation in the center of that square. 
As expected, the variation integrated over the square is smoother than the one 
in its center, with the variation being quite similar in both cases. 

In both examples shown in Figure[5] transient brightenings are visible that last 
a few minutes with an enhancement of the intensity by about a factor of about 
four. These brightenings are induced by short increases of the heating rate, which 
are ubiquit ously present in these 3D MHD models (see also iBineert and Peter , 
201lLl2013h . The brightenings shown here have properties similar to the observed 



ones in FigureEl and c o uld therefore be con s idered as a valuable model for 
blinkers ( Harrison et al . 1999; Ha rrison et al . 2003). However, the relation of 
the observed blinkers to the brightenings in these 3D MHD models will be left 
for a future study. 

In contrast to the spectroscopic observations, in the model one can compare 
the intensity variation to the magnetic field -Btr in the source region of the 
transition region emission (bottom row of Figure[5]) . It is clear that the magnetic 
field strength in the source region changes only slightly, while the intensity 
changes dramatically. The fact that t he plasma-/? parameter is much smaller than 
unity in the source region of C iv (cf . IPeter. Gudiksen. and Nordlundl . 120061 their 
Figure 12) allows very strong and rapid intensity variation while the magnetic 
field remains mostly unaffected. This opens the possibility to infer the transition 
region magnetic field despite the high temporal variability of the emission (cf. 
Section[51). 

In example (a) in the left column the magnetic field strength changes only by 
some 20%, while the intensity increases by a factor of three. With the increased 
heating the transition region moves downward to higher densities to be able 
to radiate the increased amount of heat input. Because the scale height in the 
chromosphere is only some 300 km, moving down by about 150 km is sufficient 
to relocate the transition region to densities a factor of about 1.7 higher. The 
emission goes with the density squared and thus this change of height of the 
source region by 150 km will change the intensity (and thus the radiative losses) 
by a factor of three. Over this small change in height the magnetic field is 
changing only slightly, hence the small decrease of B by some 30 G at the time 
of the peak of emission in the example (a). 

The example (b) in the right column of Figure[B]also shows only a small change 
in magnetic field while the intensity changes significantly. Here the magnetic 
field in the area is dropping towards the end of the time frame shown because of 
the (comparably slow) changes induced by the footpoint motions. This leads to 
weaker currents, less heating and consequently less transition region emission. 

In Figure|5] only J Ii os dA and the magnetic field in the source region, -Btr are 
plotted. Because according to Equation (J2|) Stokes V is basically proportional to 
/ and -Btr, and because -Btr is quite constant, Stokes V is changing in a manner 
very similar to /. This is why J|V[ os | dA is not plotted in Figure^ Normalizing 
/|Mos| dA by J I\ os dA, i.e. V as defined in Equation (0}, basically gives the mag- 
netic field, see also Section[7] and Equation ([S]). Therefore the variation of the 
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Figure 7. Spatial variation of the synthesized transition region emission, magnetic field and 
line formation height. The spatial variation is plotted along the long blue horizontal bar marked 
in in Figure^ and is shown for the same time as the snapshot displayed there. See Section r5.3l 



(normalized) Stokes V signal, V, closely follows the magnetic field in the source 
region, -Btr- Thus we plotted only -Btr but not V in Figure[6] 

The spatial variation basically shows the same properties as the temporal 
variation: strong intensity variation with only small changes of the magnetic 
field in the source region. To illustrate this, Figure[7] shows the variation along 
the blue bar in Figure[T] Despite the large-scale variation of the magnetic field 
across the region of the pore, in the region of strongest magnetic field (x from 
about 40 Mm to 50 Mm) there is only a small variation in magnetic field, while 
the C iv intensity changes by a factor of 4 to 5. Here the close relation between 
the decreased emission and magnetic field near x ~ 46 Mm to the increase in 
formation height can be seen. The small increase in formation height (due to 
less heating) by some 300 km leads to only a small change in magnetic field. Be- 
cause the (chromospheric) scale height is comparable to this change in formation 
height, the change in emission is quite dramatic. 
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Figure 8. Realistic synthetic Stokes V observation. The two panels show the synthesized 
Stokes V profiles in the 1.6 X 1.6 Mm 2 squares labeled in Figure^ accordingly. The synthetic 
observation was integrated in space over the square and in time for 12 min. The bars show a 
potential observation at a spectral sampling comparable to that of the SUMER instrument 
with a noise level of 10 -3 . See Section|6] 



6. A realistic synthetic observation: can we afford long exposure 
times? 

The above discussion shows that the observations and the synthesized emission 
from the 3D model show very strong variations in the emissivity. A future 
spectro-polarimetric instrument would have only a limited resolution in space 
and time because of the limitations in count rate and signal-to-noise to detect 
the weak Stokes V signal. To get a realistic estimate for a possible spectro- 
polarimetric observation of C iv, the estimate for the exposure time from 
Section fOl and Figure[4] adapted for specifications similar to SUMER are used. 
In particular, a resolution element of about 2" and an exposure time of 12 
minutes will be adopted. This can be considered as a worst case, because more 
modern instruments (will) have a significantly higher efficiency. If a Stokes V 
signal is visible with such a long exposure, clearly it will be visible for the shorter 
exposure times of more modern instruments. 

Figure [5] shows synthesized Stokes V profiles integrated in space over a square 
of 1.6x1. 6 Mm 2 (boxes in Figure[T]) and in time over the 12 minutes for which 
the temporal variation of the intensity is shown in these squares in Figure[6l 
To get a more realistic representation, the Stokes V spectra are shown as bars 
with a spectral sampling comparable to the SUMER instrument and with the 
addition of a noise level of 10~ 3 . 

This shows that in the both regions selected here, close to the pore and in its 
vicinity, one can detect a Stokes V signal with an instrument having a detection 
limit in Stokes V over / of about 10~ 3 . The main result from this experiment is 
that a Stokes V signal survives even when integrating in space and time, despite 
the spatial and temporal fluctuations of the intensity on scales shorter than the 
length and time scales the spectrum is integrated during the exposure. This 
is basically because the magnetic field is comparatively constant (in time and 
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space), at least much smoother than the intensity (see Section [5~3"|) . In conclusion, 
this model shows that the question posed in the heading of the current section 
can be answered with: yes, we can afford long exposure times to get information 
on the magnetic field. 

It might be that one can carry over this conclusion to existing measurements 
of the magnetic field in the upper atmosphere by spectro-polarimetry in the 
infrared. For example the signals of forbidden and o ptically-thin lines detected 
in ac tive region coronal loops with a coronagraph fe.o. lLin. Penn. and TomczvkL 
2000h might represent the background magnetic field. The polarization signals of 



the allowed He I 10830 A triplet, which is not optically thin, ob s erved in spicules 
above the limb ( Centeno. Truiillo Bueno. and Asensio Ramos . 2010h may well 
represent the magnetic field of the spicular plasma itself, as argued by these 
authors. 



7. How reliable is an inversion of the transition region magnetic 
field? 

Even if one can detect a Stokes V signal in C IV, it is not clear to what extent 
one can invert the longitudinal component of the magnetic field in the source 
region of the transition region emission. To test this, as a first step we derive 
a very simple magnetograph equation from Equation ([2]) and compare the re- 
sulting inverted magnetic field to the magnetic field in the 3D model. Using 
the wavelength and the effective Lande factor for the ClV line, A=1548A and 
<7=1.167, and assuming that the line width is comparable to the thermal width, 
AA«0.06 A, one can rewrite Equation © as 

B mc [G] w 46 000 V . (5) 

Because in the 3D model the magnetic field B is known at each grid point, one 
can now compare the inverted magnetic field B mc to the actual magnetic field 
in the source region to test this simplest inversion procedure. 

As a first step we investigate the two realistic sample spectra shown in Fig- 
ure^ Here the Stokes V signals for regions (a) and (b) are V « 0.45% and 0.55%. 
Following Equation ([5]) this corresponds to inverted magnetic field strengths of 
B nlc pa 200 G and 250 G. These inverted values are plotted as dashed lines in 
Figurc[6]on top of the variation of the magnetic field during the exposure of these 
spectra. In these two cases the simple inversion seems to work and represents 
some average value of the (line-of-sight) magnetic field in the source region of 
the transition region emission. Comparing the actual and inverted magnetic field 
in the bottom panels of Figurc|B] we estimate the error of this procedure to be 
about 20%. 

To test if the simple inversion also works in a statistical sense, we look at the 
scatter of the ratio of the inverted field B mc to the actual (line-of-sight) magnetic 
field B in the C IV source region. This is shown in Figure^ as a function of the 
Stokes V signal. At least in those areas where the Stokes V signal is above a 
noise level of 10 -3 , there is a close correlation between the inverted and the 
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Figure 9. Relation of a simple magnetic field reconstruction to the synthesized Stokes signal 
V of Civ as defined in Equation J3). Panel (a) shows the ratio of the magnetic field strength 
Bjn e reconstructed by the magnctograph equation Equation JSJ to the magnetic field at the line 
formation height, i.e., on the currogated source surface of Civ. In contrast, panel (b) shows 
the ratio of -B mo to the magnetic field at a fixed height (z = 3.3 Mm), where the horizontally 
averaged temperature is 10 5 K. See Section[7] 

actual magnetic field: the inverted signal is within about 5% of the value from 
the 3D model. However, it has to be remembered, that this is without any noise, 
which will (in part) degrade the nice correlation. 

From Figure^ it is clear that the C iv source region is far from being close 
to a flat plane. Of course, this raises the question on the usability of the in- 
verted values of the magnetic field, which mainly reflects the field on the highly 
corrugated source surface of the transition region emission in C iv as seen in 
Figure^- Therefore in Figured we compare the inverted magnetic field B me 
to the magnetic field at a constant height in the computational box. Here we 
choose a height of z=3.3 Mm, which in this model represents the height where the 
horizontally averaged temperature reaches 10 5 K. This shows a much stronger 
scatter. Still, if one assumes that the measurement would represent the magnetic 
field at a constant height, one would be able to invert the magnetic field within 
some 20% to 30% (most values in Figure^ above a Stokes V signal of 0.1% 
have a ratio from 0.75 to 1.05). Of course, this is again without assuming noise 
of the measurement. 

Certainly, the community will have to develop more elaborate inversion pro- 
cedures to interpret the Stokes V measurements in C iv in the future. However, 
even with the simple magnctograph-type inversion presented here one can hope 
to measure the transition region magnetic field within some 30%. 

8. Conclusions 

This paper presents a forward modeling of the Stokes profiles in the extreme UV 
to investigate the diagnostic potential of emission lines for investigations of the 
magnetic field in the transition region and low corona. This study shows that 
with instruments providing a polarimetric sensitivity of Stokes V/I of about 
0.1% this task to measure the transition region magnetic field can be reached; 
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probably within some 20% to 30% using simple magnetograph-type inversions 
as a first step. 

This study for C iv employs a realistic 3D MHD model that self-consistcntly 
solves for the magnetic structure and the transition region and coronal plasma 
properties. The emission synthesized from the model is highly structured in space 
and shows a strong temporal variability, just as actual observations. Naively 
one would expect that this high level of variation would destroy any significant 
polarization signal from the magnetic field. However, because the magnetic field 
structure is rather stable and smooth, the Stokes V signal survives even when 
integrating in space and time. 

Therefore spectro-polarimctcrs operating in the extreme UV will be able to 
provide reliable diagnostics of the transition region magnetic field. During its 
first flight in summer 20 10 the SUMI rock et acquired some minutes worth of po- 



larization data in Mg II ([West et all . 120111 ). A second flight took place in summer 



2012. The Civ data of this latter flight are currently calibrated and analysed, 
and it will be of high interest to see how they compare to the synthesized model 
data shown in this study — in particular if the ficld-of-vicw of SUMI would 
include some regions comparable to the magnetic structures shown here. 

Of course, it would be desirable to have measurements of the magnetic field in 
the transition region as proposed here on a space-based observatory, and embed- 
ded into a suite of instruments that will provide information on the magnetic field 



also i n the chromosphere and the corona. The SolmeX space mission (|Peter et al 



20121 ) would have provided such comprehensive measurements of the coronal 
magnetic field in the solar upper atmosphere, and was proposed to ESA recently. 
Continuing the theoretical and instrumental efforts, an opportunity might open 
up to confront the synthetic Stokes V data from the transition region presented 
here to actual observations. 
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